Violation of the Wiedemann-Pranz Law in a Single-Electron Transistor 
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We study the influence of Coulomb interaction on the thermoelectric transport coefficients for a 
metallic single-electron transistor. By performing a perturbation expansion up to second order in 
the tunnel-barrier conductance, we include sequential and cotunneling processes as well as quan- 
tum fluctuations that renormalize the charging energy and the tunnel conductance. We find that 
Coulomb interaction leads to a strong violation of the Wiedemann-Franz law: the Lorenz ratio 
becomes gate- voltage dependent for sequential tunneling, and is increased by a factor 9/5 in the 
cotunneling regime. Finally, we suggest a measurement scheme for an experimental realization. 
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Introduction. - Electron transport in conductors is 
accompanied with the transfer of both charge and heat 
(energy). Thermoelectric transport coefficients relate the 
charge and heat current, If, and Iq, to applied voltage and 
temperature differences, AV and AT, 
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The thermal conductance k is defined by Iq = kAT for 
le ^ 0, i.e., K = K - GvTS^ where S = -AV/AT = 
Gt/Gv denotes the thermopower. For macroscopic sam- 
ples of ordinary metals, the Wiedemann-Franz law pro- 
vides a universal relation between the two conductances 
by stating that the Lorenz ratio 
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(2) 



is a constant given by the Lorenz number Lq = 
{■k"^ /2>){kB / ■ It is a consequence of Fermi-liquid the- 
ory, which is applicable when screening renders Coulomb 
interactions sufficiently weak. The Wiedemann-Franz 
law indicates that both charge and heat currents are sup- 
ported by the same underlying scattering mechanisms 
with only weak energy dependence. 

The situation is fundamentally different in mesoscopic 
systems in which level quantization and Coulomb interac- 
tion drastically affects transport. The thermopower has 
been measured in_small dots with discrete level spectrum 
[H , chaotic dots [2] , carbon nanotubes and molecules ^] , 
and dots closer to the metallic (quasi-continuous) limit 
Q , and calculated for various mesoscopic systems 0, Q ■ 
Deviations from the Wiedemann-Franz law have been 
predicted for tunneling transport through quantum dots 
for weak coupling [tI, (ij , in the Kondo regime 0, 0| , for 
open dots [10], and for granular metals 11|. 

In the present Letter we address the question of 
whether and how Coulomb interaction in a metallic 
single-electron transistor (SET) with weak tunnel cou- 
plings affects the Wiedemann-Franz law. The Coulomb 
interaction plays an important role for two reasons. First, 



the finite charging energy to add or remove an electron to 
or from the island suppresses some transport processes. 
This dramatically affects the charge and thermal con- 
ductance individually, but leaves Wiedemann-Franz law 
untouched since the same transport processes are sup- 
pressed for both electric and thermal conductance. Sec- 
ond, however, Coulomb interaction leads to a strong en- 
ergy dependence of the scattering processes. This yields, 
in general, a violation of the Wiedemann-Franz law. 

By performing a systematic perturbation expansion 
of all thermoelectric coefficients up to second order in 
the tunneling conductances, we calculate the effect of 
Coulomb interaction on the Lorenz ratio. We find that 
the latter is increased due to Coulomb interaction. In 
the low-temperature regime, the Lorenz ratio becomes 
gate- voltage dependent. For the sequential-tunneling 
contributions it remains only exactly on the reso- 
nance points where the charging-energy gap of the rel- 
evant transport process vanishes. With increased detun- 
ing the Lorenz ratio rises quadratically as a function of 
the charging-energy gap. Interestingly, we find that in 
the cotunneling regime, the Lorenz ratio becomes uni- 
versal again, but by a factor 9/5 larger than Lq. 

Model. - We model the metallic single-electron tran- 
sistor by the Hamiltonian H = HL+Hn+Hi+Hch + HT- 
Here, Hr = J2kn'^kn"'lkn"'rkn with r = L,R and Hj = 
J2qn ^qnc\nCqn dcscribc noninteracting electrons in the 
left and right lead and on the island, respectively [l^ . 
Coulomb interaction on the island is described by the 
capacitance model i/ch — Ec{N — rixf' , where Ec = 
e^/(2C) defines the charging-energy scale with total is- 
land capacitance G = GL + Cn + Cg, N is the number op- 
erator of excess charges on the island, and eux = GlVl + 
GrVr + GgVg is the "external charge" that is tunable by 
gate and bias voltage. To increase the number of elec- 
trons on the island from to + 1 one has to overcome 
the charging-energy gap A^v = {N + l\Hch\N + 1) — 
{N\Hch\N) = Ec[l + 2{N - n^)]. The resonance con- 
dition Ajv = 0, where the charging-energy gap vanishes, 
is fulfilled at half- integer values of Ux- Charge trans- 
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fer processes are described by the tunneling Hamiltonian 



kqn 



-up 



h.c, where we can 



assume the tunnel matrix elements T*" to be independent 
of the states k and q and channel index n, as they vary 
on the scale of the Fermi energy, which is much larger 
than all other relevant energy scales. 

The tunnel-coupling strength for barrier r — L,R is 
characterized by the dimensionless conductance ajj — 
E„ A^r-(0)A^/(0)|T'^|2 < 1, where Njfr{0) are the den- 
sity of states of the island/leads at the Fermi level, and 



«o = J2r= 



The operator e^*'^ shifts the charge on 
the island by ±e. In general, the electron temperatures 
of left lead, island, and right lead can all be different from 
each other and differ from the lattice temperature. 

Thermoelectric coefficients. - Charge current /J — 
e{Y.kndNrkn/dt) and heat current = -{Y^kni^kr ~ 
fir)dNrkn/dt) leaving lead r — L, R are given by 

/: ^ - J [ dc.K+Hc>M + arHc<H] (3) 



= - / du;K+MC>M+a,^-HC<M] . (4) 



They depend on the rate functions for charge and heat 
transport, respectively, 



{Lo)=al / dEft{E + Lo)r{E) 



(5) 



afH=a5 / dE{E + u~ ^lr)f^{E + u;)r{E),{&) 

as well as on the correlation functions C^(w) for the is- 
land charge. The Fermi function / is denoted by /+, 
while /~ = 1 ^ /• Applied temperature or voltage differ- 
ences, AT = — Tfl and = Vl — Vr, are accounted 
for by evaluating f^{E + uj) at temperature and volt- 
age Vr, while f^{E) is taken at the island electron tem- 
perature T. Finally, we define = X]r'^e^('^)- 

In the following we concentrate on the linear-response 
regime. Furthermore, we assume that the heat current is 
conserved, i.e., the heat current entering the island from 
one lead leaves the island to the other lead and the heat 



flux from the island electrons to lattice and substrate can 
be neglected. It is convenient to use current conservation 
K/q = to write the current as h/q = {oL^Ie/q ~ 

'^o^e/q)/('^o ~^ ckQ^)' ^^'^ then expand up to linear order 
in AV or AT. Then, only the equilibrium correlation 
functions C<{uj), taken at AV = and AT = 0, enter, 
which are related to the spectral density A{uj) for charge 
excitations on the island by C^(a;) = —2TTi[l — f{uj)]A{Lu) 
and C<{uj) = 2Trif{uj)A{Lu). 

We introduce dimensionless thermoelectric 
cients gv = Gy/Gas, gr = -{e/kB)GT /G^s 
-{e/kBT)M/G^, and k = {e" /klT)K/G,,, 
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coeffi- 
m — 
where 
classi- 



cal charge conductance asymptotically reached in the 
high-temperature limit. The Lorenz ratio is L = 
{kslef [k/gv - {grlgvf] ■ We, then, find 
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with the equilibrium spectral density A{uj). 

Perturbation expansion. - We proceed by perform- 
ing a systematic perturbation expansion of the spectral 
density in the tunnel conductance ao, as was already 
done in Refs. @, based on a diagrammatic real-time 
technique [3], to address charge conductance and ther- 
mopower, respectively. To write down the first two terms 
of the perturbation expansion we only have to specify 
the results for the correlation functions from Ref. [l3| for 
vanishing voltage and temperature bias. 

By sequential tunneling processes only charge excita- 
tions at the resonances, w = A^r, can be accessed. The 
zeroth-order contribution, hence, reads 



^(0) [uj) - ^(Pjv + Pn+i)5{u - An) 



(8) 



N 



where Pjsi are the equilibrium probabilities to find 
charges on the island. The next-order contribution is 
given by A^^\uj) = X^Li ^fh"^) with 



J 



A^P{ijj) = ^^[PNa{ijj) + PN-ia+{AM + An+i- uo) + PM+ia~{AN + An-i-uj)] Rei?Ar(cj)2 
= ^^(PAr + Pw+i)<5(cc;- Aw) 



9(20Ar + (l)N+l + 4>N-l) - 



i>N+l — 



Ec 
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) + Pl 



N+1 \(PN+1 



- 4>n) 



Pn + Pn+1 



Ai\u:) = -Y,^{Pn + Pn+i)5'{u:-An)[2c^n-4>n+i~4>n-i] 



(9) 

(10) 
(11) 



Here we used the abbreviations -Rat (w) = Ajv+jO+) — AAr_i+iO+) as well as 0Ar = aoAjvRe \I'(i/3AAr/27r), 
where "^{x) is the digamma function, and d(j)N is a short no tation for dcfiM / dAjs; . 



The first term describes two-electron cotunneling that describes the charge excitations away from resonance. 
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w ^ A^r . The second and third term are contributions on 
resonance, lo = A^v, that, in the low-temperature regime, 
are identified with corrections to sequential tunneling due 
to renormalization of the tunnel-coupling strength and 
the charging-energy gap, respectively. 

Results. - Analytical results for the thermoelectric co- 
efficients can be found from Eq. © with Eqs. ((8))- pT|) . 
In Fig. [1] we show the resulting dimensionless thermoelec- 
tric coefficients gv^ Qt = itt-, and the Lorenz ratio L 
normalized by Lq, as a function of gate voltage for vari- 
ous temperatures; in Fig. [2] the temperature dependence 
of L/Lq for various gate voltages. The tunnel coupling 



is chosen as a, 



L/R 



0.01. Temperature and gate- voltage 



dependence of the Lorenz ratio can be elucidated by de- 
riving analytical expressions for various limits: 

(i) In the high-temperature regime, (3Ec ^ 1, 
Coulomb oscillations are washed out, i.e., there is no 
gate-voltage dependence. To calculate corrections to 
Wiedemann-Franz law in this regime, we expand the 
gate-voltage average of all thermoelectric coefficients in 
powers up to (/3i?c)^ to find 



L 



1 



(12) 



Deviations from Wiedemann-Franz law are visible before 
Coulomb oscillations set in, see Fig. [21 where for ksT > 
Ec the curves coincide for all gate voltages while L > Lq. 

(ii) In the on-resonance low-temperature regime, 
f3Ec > 1 but /SAat < 1 for one N (say iV = 0), trans- 
port is dominated by sequential tunneling and only the 
charge states and 1 occur. The sequential-tunneling 



contribution then yields 



L 



-^ = l + (/3Ao)V(27r)^ 

^0 



(13) 



in agreement with Ref. [8| . The Wiedemann-Franz law is 
only fulfilled for a vanishing charging-energy gap, Aq ~ 0, 
with corrections quadratic in /3Ao away from resonance. 
These corrections indicate, that the contribution of each 
transported particle to the heat current scales with the 
charging-energy gap Aq instead of temperature fc^T as 
in bulk transport. 

Higher-order corrections in ao lead to an increase of 
the Lorenz ratio. For Aq = we find 



L_ _ 1 + 4ao/3 - 2ao [7 + InjpEc/TT)] 
Lo^ 1 - 2ao [7 + HPEc/t^)] 



(14) 



with Euler's constant 7 = 0.577 . . .. The terms logarith- 
mic in temperature are associated with the renormaliza- 
tion of the tunnel- coupling strength, that enter both gv 
and k in the same way, so that it affects the Lorenz ratio 
only weakly. 

(iii) In the off-resonance low-temperature regime, 
PAff ^ 1 for all N, transport is dominated by cotun- 
neling. Expanding the thermoelectric coefficients up to 
quadratic order in temperature, and taking into account 
only Aq and A_i as the two lowest excitation energies, 
we find gv = ao(87rV3)£;^/(/3AoA_i)2, g^ ^ m ^ 0, 
and k = Sn'^/bgv- Because of weak energy dependence of 
the cotunneling scattering rate, proportionality between 
charge and heat conductance is recovered, however, with 
a different prefactor, as 



L 

To 



(15) 




FIG. 1: Coulomb oscillations of thermoelectric coefficients 
and Lorenz ratio. For high temperatures (kBT = Ec/2 - 
dotted line) oscillations are washed out. In the sequential- 
tunneling regime (ksT = Ec/W - dashed line) the Lorenz 
ratio is given by Eq. I|13[l around each resonance. For low 
temperatures {ksT = _Bc/40 - solid line) the new universal 
Lorenz ratio 9/5I/o is reached in the cotunneling regime. 
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FIG. 2: Temperature dependence of Lorenz ratio for differ- 
ent gate voltages. Two maxima separate different tunneling 
regimes. The rise to the maxima starts at ksT ~ Ec and 
ksT ~ Ao, respectively. For = 4=> Aq = Ec the two 
maxima coincide, whereas at resonance (rij, = 0.5 Ao = 0) 
the lower maximum is removed to the left. 
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FIG. 3: A possible measurement scheme. In (a) the SET 
transistor is depicted in a configuration where the drain is 
thermalized at the bath temperature and the source is a small 
conductor whose temperature can be varied. In (b) we show 
a simple thermal model of the device. 



Possible experimental realization. - A possible exper- 
imental configuration to measure the predicted effects 
is shown in Fig. [3] (a). The whole structure can be 
fabricated by standard nano-lithography and thin film 
metal deposition and with oxidized tunnel barriers in 
the SET. The normal metal source of the SET is con- 
nected to superconducting leads, some of which are 
tunnel coupled normal metal-insulator-superconductor 
(NIS) junctions and the rest are direct normal-metal- 
superconductor junctions (NS). Both the NIS junctions 
and the NS Andreev mirrors provide thermal isolation 
of the source [l^. Another role of the direct NS con- 
tacts is to suppress the charging energy of the source 
electrode while keepiiig it as small as possible to avoid 
couphng to phonons [l6|. The NIS junctions are used 
for sensitive thermometry and for heating the source; in 
proper bias range they can also cool it 15|. In contrast 
to source electrode and SET island the drain is made 
large (wide and thick) to secure thermalization at the 
bath temperature. The NS contacts are used for char- 
acterizing the SET electrically by passing current and 
measuring voltage across. The structure appears rather 
simple; we note, however, that combining a normal metal 
SET with superconducting probes requires use of non- 
conventional material combinations. Figure [3] (b) shows 
the relevant thermal model of the proposed setup. When 
applying either positive or negative heat flux into the 
source by a voltage on the NIS contacts, the dominant 
energy relaxation mechanism should be that discussed in 
this letter, i.e., electronic thermal conductance n through 
the SET transistor. This requires Ko-p k for a small 
temperature bias of the source with respect to the drain 
at the bath temperature Tq. On the left hand side of the 
inequality Kc-p = SSVTq is the linearized thermal con- 
ductance by electron-phonon coupling from source to the 
lattice. Here E 
constant 



1 • 10 WK^^m^'^ is a material specific 
15| and V is the volume of the source, which 
can realistically be made as small as 1 • 10~^^m'^ with a 
small number of probes attached. A similar argument 
applies to the choice of the size of the SET island. The 
inequality between the two heat conductances is satisfied 



at a realistic temperature of Tq = 30 mK, in particular 
near the resonance; deep in the cotunneling regime one 
may need to resort to the very different temperature and 
gate- voltage dependencies to distinguish contributions of 
K from that of Ke-p. 

Conclusions. - We have theoretically investigated 
the influence of Coulomb interaction on thermoelectric 
transport coefficients for a metallic single-electron tran- 
sistor and proposed a measurement scheme for exper- 
imental verification. We found strong violation of the 
Wiedemann-Franz law. For sequential tunneling the 
Wiedemann-Franz ratio depends quadratically on the 
charging-energy gap: the Wiedemann-Franz law is ful- 
filled only at the resonances, where the charging-energy 
gap vanishes. In the cotunneling regime the Lorenz ratio 
takes a new universal value of 9/5 of the Lorenz number. 

We acknowledge discussions with T. Heikkila and 
thank for the hospitality of the Lewiner Institute for 
Theoretical Physics at the Technion, Haifa, Israel where 
this work was initiated. 
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